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Abstract 

The model dependence of the predictions of nucleon resonances with hidden charm is investigated. 
We consider several coupled-channel models which are derived from relativistic quantum field 
theory by using (1) a unitary transformation method, and (2) the three-dimensional reductions 
of Bethe-Salpeter Equation. With the same vector meson exchange mechanism, we find that all 
models give very narrow molecular-like nucleon resonances with hidden charm in the mass range 
of 4.3 GeV < Mr < 4.5 GeV, in consistent with the previous predictions. 



1 



I. INTRODUCTION 

In the classical quark models, each baryon is made of three constituent quarks [1]. The 

pattern of the spectra and the static properties of the ground and low-lying excited states of 

baryons can be described reasonably well within these models. However, there are large devi- 

n 

ations between the predictions from these models and the experimental data [2j, such as the 
strong coupling of AT* (1535) to the strangeness, and the mass order between A/"* (1535) and 
A*(1405). In the classical 3-quark models, the iV"*(1535) with (wiid)-quarks is expected to be 
lighter than A* (1405) with (-U(is)-quarks. This problem may be solved by the penta-quark 
picture for these excited baryons. In the penta-quark models, the iV*(1535) with [mi][ds]s 
is naturally heavier than the A* (1405) with [iwZ][sg]g Actually, the conventional or- 
bital excitation energy of a original constituent quark in a baryon is already comparable to 
drag out a qq pair from the gluon field. As a result, some excited baryons are proposed to 
be meson-baryon dynamically generated states j^JloT] or states with large (qqqqq) compo- 



nents 
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121 ] . But because of the same resonances predicted by different models are in the 
similar energy region, there are always some adjustable ingredients in each model to fit the 
experimental data. Thus it is difficult to pin down the nature of these baryon resonances. 
One way to avoid such difficulty is to replace light flavor qq in these baryons by cc. Brodsky 



et al. 



131 ] proposed in the early 1980s that there are about 1% uudcc components in the 



proton. Recently, Refs.|l4Hl6|| have used different methods to predict some narrow hidden 
charm N* 5 and A* 5 resonances with masses above 4 GeV and widths smaller than 100 MeV. 
These resonances, if observed, absolutely cannot be ascribed to the conventional 3-quark 
states. Therefore, it is important to investigate the extent to which the predicted N* z and 
A* 5 resonances can be further firmly established. 

In this work we focus on the predictions [14j from a Beijing- Valencia collaboration. Their 
results are from solving the following algebraic coupled-channel equations 

T a , p (s) = Va, 7 (s)Gi(s)T^(s) + V a>f) (s) (1) 

7 

where a,/3, 7 = -DS C , DA c ,r] c N, and s is the square of the CM. energy. In Eq.© the 



meson-baryon potential is based on the vector meson-exchange mechanisms of Ref. 



5| and 
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is written in the following separable form 

\r ( \ n m P EMa + Em ? i u a\ fo\ 

V a ,p{s) = -C a , p — — -r (" ^ 0) (2) 

where the EM a is the on-shell energy of the a channel's meson, and my is the mass of 
exchange vector. The factorized propagator G 1 (E) is calculated from either using the di- 
mensional regularization or introducing a cutoff parameter A 



G(E) — > G DR {E) = f^ {a , + ln4 + m "-^ + V 4L 

167T 2 1 fx 2 2s m z B 



+ —p [ln(s - (m| - m^-) + 2gA/s) + ln(s + (m| - rrip) + 2qy/s) 
V s 

- ln(-s - (m| - m^) + 2^^) - ln(-s + (m% - m 2 M ) + 2g>/s)] } (3) 
6(B) G C{ E) = t + -»> - . (4) 

where q is on shell three momentum of MB system, u>m — \/ Q 2 + and = a/*? 2 + m^. 
It was found that the solutions of the above equations yield few narrow resonances above 
4.0 GeV. 

n 

As discussed in Ref.[14|, Eqs.([I|)-(j4]) are derived from making approximations on the 
Bethe-Salpeter (BS) equation. Schematically, the BS equation is (omitting the channel 
indices) 

T(q',q,P) = V(q',q,P) 

+ / * kV «- *• p) w^w^m^w^^ T(K p)(5) 

where P is the total four momentum of the system, q, q' and k are the relative momenta. 
For the considered vector meson-exchange, the interaction kernel is 

n*M = V*)»-< («) 

where C is a coupling constant. The complications in solving Eq.flSJ) is well known, as 
discussed in, for example, Ref.|l7| for ttN scattering. Thus approximations, such as those 
used [14j in obtaining Eqs.(lTl)- ([l"l) , are needed for practical calculations. There exist other 
approximations to solve BS equations and alternative approaches to derive practical hadron 
reaction models from relativistic quantum field theory. It is thus necessary to investigate the 



extent to which the results from Ref . [l^] depend on the approximations employed. This is the 
objective of this work. We will consider several coupled-channel models derived from using 



a unitary transformation method [18( and the three-dimensional reductions [19 1 o f Bethe 
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2l|, 



Salpeter equation. These formulations have been used in studying ttN scattering 
NN scattering [22] , and coupled-channel ttN and 7A reactions in the nucleon resonance 
region 
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24|. 



In section II, we present the considered coupled-channel formulations and discuss their 
differences with Eqs.(JT])-(|4]). The numerical procedures for solving the considered coupled- 
channel equations are described in section III. We then investigate in section IV, the numer- 
ical consequences of the differences between different coupled-channel models in predicting 
the resonance positions and the reaction cross sections. A summary is given in section V. 



II. FORMALISM 

Following Ref. [ijj , we assume that the interactions between the considered meson-baryon 
(MB) channels are due to the vector meson-exchange mechanism and can be calculated from 
the following interaction Lagrangian 

tint = £-VVV + C-PPV + £-BBV (7) 

with 

Cppv = -ig(V>[P,drP)) 

Cbbv = g((B^[V»,B]) + (B 7m B)<V» (8) 

where P and V stand for the Pseudoscalar and Vector mesons of the 16-plet of SU(4), 
respectively, and B stands for the baryon. The coupling constant g = My/2f is taken from 
the hidden gauge model with / = 93 MeV being the pion decay constant and My = 770 
MeV the mass of the light vector meson. 

By using Eq.flTj), the invariant amplitude of the PB — > PB and VB — > VB transitions 
due to the one-vector-meson-exchange interaction, as illustrated as illustrated in Fig{TJ can 
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FIG. 1: One vector exchange mechanism of meson-baryon interactions. 



be written as (suppressing the spin quantum numbers) 



M 



M 



PB,I,V 



VB,I,V 



c 



l-J 



Pv ~ m v 



C 



4/ 2 



UBilnipMi +PM j )vU B 



(9) 



<■■:.] 



4/ 2 p 2 v - ml 



UBfiAVMi +PM j )vU B j{-e* M . ■ e Mj ), (10) 



where the sub-indices i, j stand for the Mj5j and MjBj channels, I is the total isospin of the 
system, V denotes the exchanged vector meson, qi is the relative momentum of the MiBi 
channel in the center of mass frame, p a is the the four- momentum of particle a, ub 1 is 
the Dirac spinor of the baryon B^ and e^. is the polarization vector of the external vector 



meson Mi. The coefficients Cff' iy and C^ Bjy in Eqs.® and (HDD 

are taken from Ref. 



and listed in the Tables fl] and |TTJ 



TABLE I: Coefficients C^f^ B in the Eq.® for the PB system in the sector I = 1/2, 3/2, 5 = 0. 
The exchanged vector mesons V = p,oj,D* are indicated next to the values of the coefficients. 
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TABLE II: Coefficients C Mi ^^ M . B . in the Eq. (fT0|) for the VB system in the sector / = 1/2,3/2, 
S = 0. The exchanged vector mesons V = p,u,D* are indicated next to the values of the coeffi- 
cients. 







1=3/2 






P+OJ 



1=1 /2 


D*E C 


D*A+ 


J/^N 


pN 


D*Z C 


-2p+u> 





-y/3/2D* 


-1/2D* 


D*K+ 







y/3/2D* 


-3/2D* 


J/^N 


~x/zj2D* 










pN 


-1/2D* 


-3/2D* 





-2p 



We consider the coupled-channel models derived by using the unitary transformation 



method of Ref. Jjl 



181 ] and the three-dimensional reductions of Bethe-Salpeter equations 
employed in Ref. 21]. In the center of mass (CM) frame, the scattering equations within 
these models can be cast into the following general form (suppressing the spin quantum 
numbers) 



x/s - E Mk (q k ) - E Bk (q k ) + ie 



(11) 



where a = PB, VB, q, t is the relative three-momentum in channel i, y/s is the total energy, 
and E n Aai) = 



^m 2 a . + q- 2 is the energy of the particle a = M,B with a mass m ai . All 
external particles in the MB channels are on their mass-shell. Explicitly, we choose 



Qi = PMi = -PBi , 
p Mt = (E M .(p Mi ),p Mi ), 

PBi = {EBiipBi^PBi) 



(12) 



In Eq. lflTI) . the calculation of the driving term V a ' (q^, qj, y/s) from the invariant ampli- 
tude M. a,I,v (qi, qj) of Eqs.fl9|)-Eq. fri0l) and the function N(q k , y/s) in the propagator depend 



on the approximations used in deriving the above three-dimensional equations from the rel- 
ativistic quantum field theory. In the following, we specify these two ingredients in each 
model. 



A. Model based on unitary transformation method 

1. N(q k , v/i) = 1 

2. The driving term is calculated from the invariant amplitudes M. of Eqs.(l9|)-( IT0|) by ex- 
pressing the four-momentum of the exchanged vector meson in terms of the incloming 
and outgoing momenta. Explicitly, we use Eq.flH]) for PB — » PB to get 

= C M JpUB^nipMi +PM J )u\u Bl 

(PM l -PM ] r(PM l -PM :j r _ ^ {PB J -VB i r(PB r PB l Y _ 

4 u ~ v v ^ + (v "\ ^ -^ g ->- (13) 

where the kinematic variables are given in Eq. (|T2|) . 

The procedure for calculating the VB — > VB from Eq. (jT0l) is same since the factor 
—e* M . ■ eyij does not depend on p v of the exchanged vector meson. 

B. Models based on three-dimensional reductions 

We consider the three-dimensional reductions developed by Kadyshevsky, Blankenbecler 



and Sugar, and Thompson, as explained in Ref . 21] . All have the same form of the potential 
which is defined by setting the time component of the four-momentum of the exchanged 
vector meson V to zero. By using Eq.(j9]) for PB — > PB, we get 

,PB,I,vM v U Bj [ Ml {EmAPmJ + E Mi {p Mj ))r - iPMi + PMj) ■ l\u Bi 



c 



4/ 2 -pi-mi 

(14) 



The procedure for calculating the VB — > VB from Eq.f fTOl) is same since the factor — e* M . -cmj 
doe not depend on p v of the exchanged vector meson. 

The function N(k, y/s) of the propagator of Eq. ffTTl) for each reduction is 



1. Kadyshevsky: 



N( qi ,y/s) = l 



(15) 



2. Blankenbecler-Sugar : 



3. Thompson: 



at( ^ 2(E Mi ( qi ) + E Bi ( qi )) 

yfs + {E M .( qi ) + E B .( qi )) 



N{ qi ^) = ^M±^m. (it) 
v s 



Note that for the on-shell momentum, defined by y/s = E M .( qoi ) + E B .( qoi ), N( qoi , Vs) = 
1. Thus all models satisfy the same unitarity condition defined by the cuts of the propagators 
of Eq.flUJ. 



III. CALCULATION PROCEDURES 

In this section, we describe our procedures for solving the coupled-channel equations to 
obtain the MB — > MB cross sections. 

It is convenient to cast Eq.f fTTj) into the following familiar form 

+ Y] / d qk V a '\ quqk ,^)^ — r-—^——T a '\ qk , qj ,^Ts) 

(18) 

where a = PB, VB, and 

V^( qi , qj , V~s) = iV 1/2 (%, Vs) V^m, qj , v^)iV 1/2 (g„ Vs) (19) 
T^( qi , qj , Vs) = N 1 / 2 ^, V^)f aJ m, qj , VslN 1 ^, V~s) (20) 

With the normalization < p\p' >= S(p—p) for plane wave states, we obtain from 
Eq.( fl8|) the following coupled-channel equations in each partial wave 

T iis x m A (QuQ2,Vs) = V£* Sl M , S2 ( q i, q2 ,Vs) 

+ J Qs^V&uL^sM' ?3, Vs) G( q3 , Vs) T^ S3jL2jS2 ( q3 , q 2 , Vs) 

(21) 



where J is the total angular momentum; Li and Si are the orbital angular momentum and 
total spin of the MiBi channel, and the propagator is 



The matrix elements of the potential in E g .([2 ID can be conveniently calculated from Eg. (119 
by using the LSJ-helicity transformation [23[. Explicitly, we obtain 

V /(2L 1 + 1)(2L 2 + 1) 1 / m Bl m Bi 

— - t L 1 ,L 2 {Ql,Q2) 



2 J + 1 (2tt) 3 V 2E M (q 1 )E B (q 1 )2E M (q 2 )E B (q 2 ) 

V^^y/.V \^ \ ^ r JM Sl r Si,M Sl r JMs 2 r s 2 ,Ms 2 

X °"1,2 O ii,5i,0,A/s 1 SAf 1 AM 1 Js 1 -AB 1 O L2,5 2 ,0,Af S2 O jAf 2 AM 2 ,iB 2 -AB 2 
V Ajv/jAs-j^ Am 2 Ab 2 

xiV 1 / 2 ^; v^) < gi;-A Bl ,A Ml |V J |A M2 ,-A B2 ;g 2 > iV 1/2 (g 2 , y/7s) (23) 

where A Q is the helicity of particle a, and by writting Eq. flTB"]) or Eq. (fl4|) (and also the 
similar forms for VB) in helicity representation we can evaluate 

< Qi', — Abd ^Mi\V j \Xm 2 j ~ As 2 ; q 2 > 



(24) 



where cos# = g\ • q 2 , and the matrix element in the integrand can be calculated by writing 
Eq.( ll3p or Eq. (1141) (and also the similar forms for VB) in helicity representation for each 
of the considered coupled- channel models. 

In Eq. (j23~j) Cj^. j 2>mj2 =< J, M\ji, j 2 ,rrij 1 ,mj 2 > is the Clebsh-Gordon coefficient, the 
isospin factor is 



m 'M 1 m I Bl m lM 2 m lB 2 



where (ImTtiim, lBTni B ) are the isospin quantum numbers of MB, and the form factor is 
chosen as 



At/ n\, At 



= (t^^^t^)^, (26) 

where the cut-off parameter Ay is assumed the same value for all exchanged vector mesons 
for simplicity. 



The differential cross sections are calculated from the partial-wave amplitudes by 



da = E Ml E B E, h E B ^ ^ l<MlB2lT(V?)lMlBl>n27) 

dil s qi (2j Ml + l){2j Bl + 1) ^ ^ 



with 

<M 2 B 2 \ T (y/s)\M x B x > 

J, I L\,Si,L2,S2 

U Li,Si,0,Msi3Mi >3Bi,mj Mi ,m Jfli ^ L 2) S 2 ,Ml 2 Ms 2 3m 2 ,3b 2 > m 3 M2 ' m 3B 2 
X C hi 1 ,Ibi ,m lMi ,m lBi C I M2 Jb 2 ,mr Ma ,m lB ^ ■ ( 28 ) 

Obviously, Mj = M Sl = m jMi +m jBi , Mg 2 = m jM2 +m jB2 , M L2 = (m jMi +m jBi ) - (m jM2 + 
m jB ) and Mj = m I& , + m /si = m H , + m Ig . 



IV. THE RESULTS 



In this section, we show the results from four models listed in SecJTTl and then discuss 
their differences with previous works [1 



A. The results of 4 models listed in SecJlTl 

We first consider the model based on the unitary transformation method described in 
subsection II. A. We determined the resonance pole positions (Mr = M — i^) by using 
the analytic continuation method of Ref. [25[. We find that the resonance positions are 
sensitive to the cutoff A, as seen in Table IIHI The resonances are generated only when the 
cutoff is larger than 800 MeV. As the cutoff A changes from 800 MeV to 2000 MeV, the 
"binding energy" (AE = M - E thr ) is changed greatly from 0.002 MeV to 23.9 MeV. The 
corresponding changes in imaginary parts are also very large. 

These resonances are very close to the threshold of -DS C in the PB sector and -D*S C 
in the VB sector. They are mainly caused by the strong attractive potential from the t- 
channel p meson exchange in DH C — > DH C and D*H C — > D*S C processes. The situation 
here is different from the case when only light flavors are involved. For the ttN interaction, 

10 



there is no resonance below the irN threshold, although the t-channel p meson exchange also 
provides attractive potential there with a similar coupling constant. The differences between 
two cases are mainly from the term (pM t +Pa/ 3 ) i n Eq. (jl3j) . The potential is proportional to 
(rriMi +^Mj) near the threshold of system. For the -DE C case, (m^. +m^) ~ 4 GeV, while 
it is about 0.3 GeV for the ttN case. Hence the attractive potential of D£ c is an order of 
magnitude stronger than that of irN. This give a natural explanation why the PB system 
with heavy quarks can have quasi-bound states while the corresponding pure light quark 
sector cannot. The similar thing happens also for the VB system. 

TABLE III: The pole position (M - iY/2) and "binding energy" (AE = E thr - M) for different 
cut-off parameter A and spin-parity J p . The threshold Ef^r is 4320.79 MeV of DT, C in PB system 
and 4462.18 MeV of D*Y, C in VB system. The unit for the listed numbers is MeV. 







PB System 




VB System 




Jf = 




M-iT/2 


AE 


M-iT/2 


AE 




650 












800 






4462.178 - 0.002i 


0.002 




1200 


4318.964 - 0.362i 


1.826 


4459.513 - 0.417i 


2.667 




1500 


4314.531 - 1.448* 


6.259 


4454.088 - 1.662z 


8.092 




2000 


4301.115 - 5.835i 


19.68 


4438.277 - 7.115* 


23.90 


JP = 


3~ 
2 












650 












800 






4462.178 - 0.002i 


0.002 




1200 






4459.507 - 0.420i 


2.673 




1500 






4454.057 - 1.681* 


8.123 




2000 






4438.039 - 7.268i 


23.14 



The three other models based on three-dimensional reductions in subsection II. B give 
similar results as shown in Table IIVI together with those from the model based on unitary 
transformation method, taking the cut-off parameter A = 1500 MeV. The corresponding 
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results for the total cross section of r] c p — > r] c p are shown in FigfSJ All these four models 
predict a resonance below the -DS C threshold. The masses and widths of the resonances 
from these different models are almost the same. 

TABLE IV: Comparison for 4 models with the cut-off A = 1500 MeV and J p = 1/2" for PB 
system, where the threshold energy E t h r is 4320.79 MeV of Z)£ c . "A" is for the model based on 
unitary transformation method; "B" is for Kadyshevsky model; "C" is for Blankenbecler-Sugar 
model; "D" is for Thompson model. AEa and Ta are the binding energy and width for the case 
A. The unit is MeV. 



Models 


M - 


iT/2 


AE 


1 AE-AEa I 
1 AE A 1 


ir-r A , 
1 r A 1 


A 


4314.531 


- 1 Am 


6.259 








B 


4314.983 


- 1.737? 


5.807 


7.222% 


19.96% 


C 


4314.436 


- 1.879? 


6.354 


1.518% 


29.77% 


D 


4314.824 


- 2.04H 


6.966 


11.30% 


40.95% 




FIG. 2: The total cross section of r] c p — > r] c p vs CM. energy is shown for 4 models. The 4 lines 
correspond to 4 models listed in Table HVl 
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B. Comparison with previous works 

In Ref.jjjJ using Valencia model, the mass and width of predicted resonance in the PB 
system is about 4265 MeV and 23 MeV (for rj c N channel only). Both binding energy and 
width are much larger than the results in this work. All models considered in this work differ 
from the model used in Ref. 14 in calculating the MB — > MB potentials Eqs.(!9|)-( TT0|) . We 
will get the form used in Ref . |14 1 . if we : (1) neglect the lower component of Dirac spinor and 
keep only the time component 7 ; (2) set the momentum squared of the exchanged vector 
meson V to be p v = (E^_ — E^) 2 — (q° n — q° n ) 2 , where the E^. and q° n are, respectively, 
the on-shell energy and momentum of the meson in channel i. We have investigated the 
effects from taking each of these two assumptions. If we only make the first simplification 
by neglecting spin of baryons, then the corresponding results for the resonance are shown as 
for potential A' in Table |Vj If we continue to make the second simplification by neglecting 
the momentum of the exchanged vector meson, the results are shown as for potential A" of 
Table [V] In Figj3j we show the results corresponding to these two simplifications (dashed 
line for A' and dot-dashed line for A") for the r\ c p — > rj c p total cross section, compared with 
that (solid line for A) from the model based on unitary transformation. Clearly, the second 
simplification shifts the resonance position to a much lower value and also increases the width 
significantly. This is the main reason for the difference between our present results and those 
from Ref. 14J. The second simplification makes p v and hence the potential V independent 



on the integral momentum in Eq.(j5]) so that Eq.([5]) is simplified to Eq.([T]) instead of Eq.( fl8l) 
where the potential V with integral momentum dependence is inside the integration. Eq. ([T]) 
and Eq. ffTSl) give the different results. 
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TABLE V: Comparison for different potential approximations with cut-off A = 1500 MeV and 
J p = 1/2- for PB system. The threshold E thr is 4320.79 MeV of DT, C . "A" is for the full 
potential; "A"' is for the neglect of spin of baryons; "A"" is for the neglect of both spin of baryons 
and momentum of exchanged vector meson. AEa and Ta are the binding energy and width for 
the case A. The unit is MeV. 



Potential 


M - 


*T/2 


AE 


1 AE-AE A | 
1 AE A 1 


1 T - T A 1 

1 T A 1 


A 


4314.531 


- 1.448i 


6.259 








A' 


4316.315 


- 0.967i 


4.475 


28.50% 


33.22% 


A" 


4229.362 


- 3.914i 


91.43 


1361% 


170.3% 




FIG. 3: The total cross section of rj c p — > rj c p vs CM. energy for different potential approximations. 
The solid line is for the full potential, corresponding to A in Table |Vj The dashed line is for the 
neglect of spin of baryons, corresponding to A' in Table |Vj The dot-dashed line is for the neglect 
of both spin of baryons and momentum of exchange vector, corresponding to A" in Table |V] 



f.Q, 



Besides Ref . [14| , there are two later publications |15|, Il6] also predicting the existence of 
A^* around 4.3 GeV with hidden charm. 

In Ref. lol]. the S-wave S C D and A C D states with isospin 1=1/2 and spin S=l/2 are 
dynamically investigated within the framework of a chiral constituent quark model by solving 
a resonating group method (RGM) equation. The calculation not only includes vector 



14 



mesons (p and uS) exchange, but also scalar(cr) meson exchange, which provides an additional 



attractive force. Therefore, the binding energy in Ref. 



15] is larger than that in this work. 



The mass of the bound state of -DS C is about 4279 — 4316MeV. 



In Ref. 16], the Schrodinger Equation was used to find the bound state of -DS C and -D*S C 
with effective meson exchange potential. For the PB system, the p, u and a exchanges were 
considered. They tried the different sign of coupling constants of various vertices. When 
they chose the u exchange to be repulsive, and the p and a exchange to be attractive, they 
also found the isospin 1/2 bound state of -DE C with cut off A > 1.6 GeV. The binding energy 
is about — 16MeV corresponding to A = 1.6 — 2.2GeV, similar to the results in this work. 

V. SUMMARY 

We have investigated the possible existence of nucleon resonances with hidden charm 
within several coupled-channel models which are derived from relativistic quantum field 
theory by using a unitary transformation method and the three-dimensional reductions of 
Bethe-Salpeter Equation. With the same vector meson exchange mechanism, we find that all 
models give very narrow molecular-like nucleon resonances with hidden charm in the mass 
range of 4.3 GeV < Mr < 4.5 GeV, in consistent with the previous predictions. From our 
analysis, the heavy mass of particles with the core components would make the attractive 
potential stronger than the case with only light flavors. The widths of these resonances are 
very narrow in our models, because they need heavy vector meson D* exchange to decay 
to open channels. Furthermore, we compare our results with previous works. All of models 
predict a resonance below the DH C threshold. We also find that the pole position would be 
shift a lot if we set p v = {E°j^_ — E^.) 2 — (q° n — q° n ) 2 for the exchanged vector meson V in 
the potential. We look forward to find these predicted resonances with hidden charm in the 
reactions, such as e p —> e J/ip p, p p — >■ p r] c (J/ip) p, and p p — >■ p r] c (J/ip) p. 

Similarly the super-heavy N* with hidden beauty should also exist although the binding 



energies may be not as large as given by the simple Valencia model calculation of Ref. 



26] 
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